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ABSTRACT 


Tests  have  been  made  of  the  Runge-Kutta  and  (especially)  Nunerov 
methods  of  solving  the  reduced  Schr6dinger  wave  equation  for  a  tabulated 
Morse  potential.  This  is  preliminary  to  numerical  solution  of  the 
equation  for  a  Klein- Dunham  potential. 

Normalized  wave  functions  for  ,  B"Vl,  of  orders  0  and  5*  given  at 
an  interval  of  0.01A  in  internuclear  separation,  agree  closely”  with  those 
obtained  analytically. 

A  criterion  fer  goodness  of  wave  functions  generatod  numerically, 
depending  on  ob&orvsd  rotational  constants,  is  described. 


1 .  INTRODUCTION 


This  report  describes  the  first  phase  of  a  project  to  solve  the 
reduced  SchrSdinger  wave  equation  numerically,  given  a  tabulated  Klein- 
Dunham  potential  energy  function  (Klein  1932,  Dunham  1932,  Jarmain  I960, 
196l).  The  'true'  potential,  as  is  well  known,  is  not  available  ir 
closed  form. 

For  certain  diatomic  molecular  transitions,  Morse  vibrational  wave 
functions  are  inadequate  for  the  calculation  of  Pranck-Condon  factors 
and  other  parameters.  Therefore  we  wish  to  determine  how  large  the 
discrepancies  are  that  sometimes  exist  between  Franck-Condon  factors 
based  on  a  Morse  potential,  and  those  based  on  a  more  accurate  Klein- 
Dunham  function. 

Such  factors  are  often  useful  in  their  own  right  in  predicting 
approximate  relative  intensities  of  molecular  bands,  and  of  course  have 
been  used  extensively  in  conjunction  with  good  measured  intensities  to 
yield  the  shape  of  the  electronic  transition  moment  (see  o.g.  Hebert 
and  Nicholls  1962). 

The  present  summary  deals  only  with  the  numerical  solution  of  the 
wave  equation  for  a  tabulated  Morse  potential.  It  was  considered 
prudent  first  to  verify  that  correct  Morse  wave  functions  could  be  so 
generated,  especially  -when  these  were  available  for  checking  from 
earlier  analytic  work. 


2 .  Solution  of  the  Differential  Equation 

Neglecting  rotation,  the  reduced  Schroainger  equation  satisfied  by 

a  vibrational  wave  function  ^  (r) ,  with  v  the  vibrational  quantum  number, 
^2^  ^  —  __ 

— ^  +  SlUi  Ev  -  U(r)  I  =  0  (1) 

dr  h  ]_  j 
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where  r  is  internuclear  separation,  ji  is  the  reduced,  mass  of  the  molecule, 
E^_  is  the  appropriate  eigenvalue,  and  U(r)  is  the  vibrational  potential 
energy  function. 

Tests  of  the  method  described  below  were  carried  out  on  IF, ,  B^n  for 
the  following  reasons.  Firstly,  ordinates  of  Morse  wave  functions  up 
to  order  12  had  been  calculated  previously  for  this  state  from  formulae 
(Jarmair.  and  ITicholls  1952.  I960).  Secondly,  there  is  good  agreement 
(mostly  less  than  1$  deviation)  between  the  tabulated  Morse  and  Klein 
energies  at  specific  values  of  r.  Finally,  if  a  theoretical  rotational 


constant  a 
relation 


.•0 


3  =  B  -  a  (v  +  i) 

v  e  e  2' 


6  E 


,  is  determined  from  Pekeris'  (1934) 

(2) 


CO  X 


c  e  _j 

it  agrees  with  the  observed  a  to  about  0.5 Thus,  anticipating  the 
use  of  a  Klein- Dunham  potential  to  obtain  wave  functions,  one  would 
expect  fairly  good  point  by  point  agreement  with  corresponding  Morse 
functions.  This  would  provide  a  preliminary  check  on  the  validity  of 
the  procedure  used  to  start  the  Kloin-Dunham  solution.  A  later  report 
will  discuss  that  procedure. 

Solutions  were  first  obtained  using  a  Punge-Kutta  fourth  order 
process  as  adapted  for  digital  computers  by  Gill  ( 1 9 5 4 ) •  This  choice 
was  made  because  it  was  available  as  a  subroutine  in  the  U.W.O.  Computing 
Centre  Library,  because  it  is  known  to  be  very  stable  and  requires 
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knowledge  of  only  one  point  and  a  slope  at  that  point  for  starting. 
These  values  wore  determined  from  analytic  Morse  expressions,  on  the 
extreme  left  side  of -the  wave  function  beyond  which  point  ordinates 
are  so  small  that  they  contribute  nothing  significant,  for  four  figure 
accuracy,  to  an  overlap  integral  involving  a  wave  function  of  another 
electronic  state.  Integration  then  proceeded  from  left  to  right. 

As  anticipated,  this  method  was  effective  but  relatively  slow. 

It  did,  however,  indicate  sensitivity  of  the  curve  to  choice  of  eigen¬ 
value  and  also  a  suitablo  sub-interval  size  for  integration.  For  each 
of  tho  two  orders  tested,  a  start  was  made  with  the  experimental  eigen¬ 
value,  but  minor  adjustment  in  it  was  necessary  to  keep  ordinates 
numerically  small  over  a  reasonable  interval  on  the  right  hand  side. 

At  best,  wave  functions  remained  unbounded  at  r  =  ,  but  were  simply 

cut  off  in  the  region  of  smallest  ordinates.  hTot  surprisingly,  the 
subinterval  had  to  be  diminished  as  the  order  of  tho  wave  function 
increased . 


The  Runge-Kutta-Gill  routine  was  then  abandoned  in  favour  of  a 
fonr.ul"tion  due  to  Numerov  (  1933)  a  partial  derivation  of  which  follows. 
Taken  fror  central  difference  theory,  equation  (3)  (Hartree  "'958)  relates 
the  second  difference  of  a  function  to  its  second  derivative  and  differ¬ 
ences  of  the  latter,  namely, 


-  (irM *6'  +  \ 


Equation  (l)  to  be  solved  is  of  the  form 

4> '  '  =  F  (r)  l> 


- —  1  ' 

240  v0  i 


+  0  (6r)q  (3) 


(4) 
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Neglocting  the  fourth  difference  in  (3),  we  have 


\  “  [*0*  +  12  *2  ^c/j 


It  is  also  true  that 


*  4^  =  4»1  -  24*0  +  4'_1  (6) 

where  the  subscripts  refer  to  consecutive  values  of  Again,  by  virtue 


of  (4) 


*o'  -  F0  ^0 


'  •  _  F  (L  _  2^  1  +  F  1 

9  ^0  1  *1  0  ^0  -1  -1 


Therefore,  combining  equation0  (5).  (6)  and  (7), 

[1  -  4(V)2  f]  -  2  [l  -  -(Srf  F0]  *0  +  [1  -  ±(Sr)2  fJ*.,  -  (irf  ^ 

(8) 

Finally,  a  formula  for  4 ^  in  terms  of  Lq  and  is 

0  -  h{^)2  F]  h  •  [2  *  I  ('r)2  Fo]  *0  -  -  T2(4r)2  p-3*-i  (9) 

Thus.  Numerov’s  equation  ',9)  generates  the  functional  value  if 
two  starting  values  and  4>  ^  arc  provided,  and  of  course  the  calculated 
4^  is  needed  together  with  1q,  to  produce  at  the  next  step,  and  so  on. 

The  technique  is  fast  and  sufficiently  accurate  for  the  purpose,  providing 
a  suitably  fine  subdivision  is  used.  As  mentioned  before,  the  observed 
eigenvalue  does  not  give  a  wave  function  that  is  well  behaved  on  the 
right  hand  side.  This  is  probably  mainly  due  to  the  approximate  nature 
of  the  procedure. 


3.  RSSULTS 

Morse  wave  functions  of  order  0,  5  and  10  have  been  constructed  using 


C, 
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TA3LE  I 


MORSE  WAVE  FUNCTION 
n2  B3n  v  =  0 


r(A) 

Analytic 

Numerical 

r(A) 

Analytic 

Numerical 

1.02 

1 .011x10 

1 .011x10 

1.23 

3.149x104 

3.149x104 

1.03 

2.410x10 

2.410x10 

1 .24 
'# 

2.944xl04 

2.944x104 

1.04 

5*430x10 

5.430x10 

1.25 

2.662x104 

*2  .66lx104 

1.05 

1 .158x102 

1.1 58x1 O2 

1 .26 

2.328x104 

2.328x104 

1.06 

2.340x102 

2. 340x1 O2 

I.27 

1 »973x104 

1 -973x104 

1 .07 

4-488x102 

4-488x102 

1.28 

1 . 620x1 04 

1.62 0x1 O4 

1.08 

8.l79x102 

*8-180x102 

1.29 

1 .290x104 

1 .290xl04 

I.09 

1 .418x103 

1 .418x103 

1.30 

9-975x103 

9-975x103 

1.10 

2.342x103 

*2 . 343x1 03 

1.31 

7.491x103 

7.491x103 

1.11 

3 .690x103 

3.690x103 

1.32 

5.469x103 

5-469x103 

1.12 

5-5^9x103 

5. 549x10 3 

1.33 

3.883x103 

3.883x103 

1.13 

7-977X103 

7  -977x103 

1 . 34 

2.684x103 

2.684x103 

1.14 

1 .O97XIO4 

1 .097XIO4 

1.35 

1 .8O7XIO3 

I.8O7XIO3 

1-15 

1 .446x"04 

1  .446x1 04 

1.36 

1.186x103 

1 .186x103 

1.16 

1 =827x104 

1  .827x1 04 

1.37 

7-586x102 

7  «586x102 

1.17 

2.215x10‘t 

2  .2 15x1 04 

1.38 

4.737XIO2 

*4*736x102 

1.18 

2.582x104 

2. 582x1 04 

1-39 

2.888x102 

*2.885x102 

1519 

2. 893x1 O4 

2. 89 3x1 04 

1.40 

1  -720x102 

*1-7l6x102 

1.20 

3.122x104 

3.122x104 

1.41 

1 .002x102 

*9*941x10 

1.21 

3.245x104 

3. 245x1 04 

1.42 

5.7O4XIO 

*5.577x10 

1.22 

3.254x104 

3.254x104 

Asterisks  indicate  where  analytic  and  numerical  values  differ. 
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TABLS  II 

MORSE  HAVE  FUNCTION 
N0  B]n  v  =  5 


r(  A) 

Analytic 

Numerical 

r(  A) 

Analytic 

Numerical 

i  .01 

5*734x102 

*5-735x102 

1.29 

-1.600x103 

-1 .600xl03 

1.02 

1 .187x103 

*1.188xl03 

1.30 

8.854xl03 

*8.855xl03 

1.03 

2.285x103 

*2.286x103 

1.31 

1 .677XIO4 

1 .677x10" 

1 .04 

4. 08 3x1 03 

4. 08 3x103 

1.32 

2  . 0 1 1  x  1 04 

2 .01 IxlO4 

1.05 

6 .76OXIO3 

6 .76OXIO3 

1.33 

1  .83lxl04 

1 .83lxl04 

1.06 

1 .034xl04 

1. 034x1 O4 

1.34 

1.215x104 

1 .2 1 5x104 

1 .07 

1 . 455x1 O4 

1 .4 55x1 04 

I.85 

3.308x103 

3. 308xl03 

1.08 

1  .869x104 

1 .8  69x1 04 

1.36 

-6.216x103 

-6 .2l6xl03 

1.09 

2.l67x104 

2 . 167x10" 

1.37 

-1.468xl04 

-1. 468x1 O4 

1.10 

2.220xl0/r 

2.220x10' 

1.38 

-2 .09 IxlO4 

-2. 09 IxlO4 

1.11 

1 .921x104 

1 .92 IxlO4 

1.39 

-2.442x10 

-2.442x10" 

1.12 

1 .2  44x1  O'" 

1.2/>4x104 

1-40 

-2. 5 30x1 04 

-2. 530x1 04 

1.13 

2.802x103 

2.802x103 

8.41 

-2.406x10" 

-2. 406x1 04 

1.14 

-7.524x103 

-7.524x103 

1.42 

-2.<39x104 

-2.l39x104 

1.15 

-1.564x104 

-1 .564x1 o4 

1.43 

-1 .797XIO4 

-1 .797XIO4 

1.16 

-1 .887XIO4 

-1 . 887x10" 

1.44 

-1 .436x104 

-1.438x104 

1.1? 

-1.592xl04 

-1. 592x1 04 

1.45 

-1.102x104 

- 1 . I02xl04 

1.18 

-7.539x103 

-7.539X103 

1.46 

-8 . 1 14x103 

-8 • 1 14x1 O3 

1.19 

3»48lxl03 

3.48 IxlO3 

1.47 

-5.765x103 

-5 .765XIO3 

1.20 

1 .32 4x1 O4 

1. 324x1 O4 

1 .48 

-3.96lxl03 

-3.96lx103 

1.21 

1 .8l7x104 

1 .817X104 

1.49 

-2  .  638xl03 

-2.638x103 

1.22 

1 . 645x1 O4 

1.645x10" 

1 .50 

-1.706x103 

-1.706xl03 

1.23 

8*723x103 

8-723x103 

1-51 

-1 .073x103 

-1.073x103 

1.24 

-2.l75xl03 

*-2 . 176x103 

1-52 

-6 .577XIO2 

*-6 . 576XIO2 

I.25 

-1.23lxl04 

-1.23lxl04 

1.53 

-3.930x102 

*-3.929xl02 

1.26 

-1.8l4x104 

41.814x104 

1.54 

-2.293xl02 

*-2.290x102 

1.27 

1.28 

-1 *779x104 
-1.154xi04 

-1 J779x104 

-1.154xi04 

1.55 

-1.308x102 

*-1 .302x102 

Asterisks  indicate  where  analytic  and  numerical  values  differ. 
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an  I.B.K.  65O,  with  detailed,  results  given  for  the  first  two.  The 
tenth  order  is  still  under  study,  with  indications  that  a  sub-interval 
of  loss  than  0.005A  is  noodod.  Intervals  used  for  the  others  were 
0.01A  and  O.OO5A  respectively. 

Table  I  shows  a  coup a risen  between  normalized  ordinates  obtained 
analytically  to  four  figures  and  those  from  the  present  method,  for 
v  ■  0.  The  final  eigenvalue  was  863.44463  cm  ^  compared  with  the  ob¬ 
served  863.438  cm  '.  Starting  values  at  1.02A  and  1.03A  were  10.110 
and  24*097  respectively,  with  r  measured  in  cm.  Table  II  gives  a 
similar  comparison  for  the  fifth  order,  with  eigenvalue  of  9098*3955  cm  ^ 

against  the  observed  9099*388  cm  ' .  Starting  values  at  1 .010A  and 

2  2 

1.0l5Awere  5*7332  x  10  and  8*3259  x  10  respectively.  In  general, 
agreement  is  perfect  to  four  figures  except  on  the  extreme  right. 

4.  B  Test 

v _ 

It  was  felt  from  the  beginning  that  as  many  chocks  as  possible 
should  be  provided  for  the  normalized  Klein- Dunham  wave  functions 
eventually  to  be  constructed  numerically.  Otherwise  it  would  be  un¬ 
safe  to  claim  that  these  functions  are  necessarily  an  improvement  on 
Morso.  One  tost  that  is  immediately  available  is  of  course  orthogonality 
between  wave  functions  of  different  order.  A  second  one,  which  depends 
on  observed  rotational  constants,  was  devised  by  Dr.  P.A.  Fraser. 

According  to  Herzbcrg  (l95^),  P*  106,  it  is  plausible  to  use  a 
mean  value  for  the  rotational  constant  B  in  the  vibrational  state  con¬ 
sidered,  that  is 
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where  r  is  the  mean  value  of  — -  during  the  vibration.  See  also 

c-  £  /  n 

r  r  T 

Hyllcraas  (1936)  for  justification.  But  ■'r  is  given  as  well  by 


which  can  be  evaluated  numerically  for  any  tabulated  wave  function  over 
the  significant  range  of  r.  Note  that  the  Morse  wave  function  is  not 
zero  at  the  origin*  although  extremely  small,  and  therefore  the  integrand 
is  infinitely  great  at  the  lower  limit.  In  practice,  however,  we  in¬ 
tegrate  no  noaror  to  the  origin  then  say  O.5A  (more  in  this  case),  because 
ordinates  arc  then  trivially  small.  For  purposes  of  the  test,  one  can 
assume  over  the  remaining  interval  a  segment  of  the  true  wave  function 
which  will  be  of  a  high  enough  order  of  smallness  (zero)  at  the  origin 
to  overcome  the  effect  of  r  ^  in  the  integrand.  Thus,  comparison  of  the 
right  hand  side  of  equation  (lO)  with  the  experimental  should  give  a 

good  indication  of  the  correctness  of  the  wave  function  1/  . 

v 

Results  of  applying  this  test  to  the  Morse  functions  above  were 

—  *  —  A 

Bq  =  I.629O  cm  compared  with  1 .628°  cm  1  obser/ed,  and  3,.  =  1.5362  com¬ 
pared  with  1.5368.  The  wave  function  of  lower  order  is  evidently  more 
nearly  correct  on  the  basis  of  the  above  criterion.  It  is  hoped  that 
further  tests  can  be  developed  as  -work  continues. 


5.  FUTURE  PROGRAM 

The  next  report  will  deal  with  construction  and  testing  of  Kleirr* 
Dunham  wave  functions  obtained  by  the  present  method  for  ,  B^TT  and 
possibly  other  molecular  states.  This  work  is  already  under  way  but 
subject  to  revision  and  therefore  not  roady  for  publication. 


-  o  - 
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